Identification of repurposing therapeutics toward SARS-CoV-2 main protease by virtual screening.

SARS-CoV-2 causes the current global pandemic coronavirus disease 2019. Widely-available effective drugs could be a critical factor in halting the pandemic. The main protease (3CLpro) plays a vital role in viral replication; therefore, it is of great interest to find inhibitors for this enzyme. We applied the combination of virtual screening based on molecular docking derived from the crystal structure of the peptidomimetic inhibitors (N3, 13b, and 11a), and experimental verification revealed FDA-approved drugs that could inhibit the 3CLpro of SARS-CoV-2. Three drugs were selected using the binding energy criteria and subsequently performed the 3CLpro inhibition by enzyme-based assay. In addition, six common drugs were also chosen to study the 3CLpro inhibition. Among these compounds, lapatinib showed high efficiency of 3CLpro inhibition (IC50 value of 35 ± 1 µM and Ki of 23 ± 1 µM). The binding behavior of lapatinib against 3CLpro was elucidated by molecular dynamics simulations. This drug could well bind with 3CLpro residues in the five subsites S1', S1, S2, S3, and S4. Moreover, lapatinib's key chemical pharmacophore features toward SAR-CoV-2 3CLpro shared important HBD and HBA with potent peptidomimetic inhibitors. The rational design of lapatinib was subsequently carried out using the obtained results. Our discovery provides an effective repurposed drug and its newly designed analogs to inhibit SARS-CoV-2 3CLpro.

Computational study on peptidomimetic inhibitors against SARS-CoV-2 main protease.

The emergence outbreak caused by a novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has received significant attention on the global risks. Due to itscrucial role in viral replication, the main protease 3CLpro is an important target for drug discovery and development to combat COVID-19. In this work, the structural and dynamic behaviors as well as binding efficiency of the four peptidomimetic inhibitors (N3, 11a, 13b, and 14b) recently co-crystalized with SARS-CoV-2 3CLpro were studied and compared using all-atom molecular dynamics (MD) simulations and solvated interaction energy-based binding free energy calculations. The per-residue decomposition free energy results suggested that the key residues involved in inhibitors binding were H41, M49, L141-C145, H163-E166, P168, and Q189-T190 in the domains I and II. The van der Waals interaction yielded the main energy contribution stabilizing all the focused inhibitors. Besides, their hydrogen bond formations with F140, G143, C145, H164, E166, and Q189 residues in the substrate-binding pocket were also essential for strengthening the molecular complexation. The predicted binding affinity of the four peptidomimetic inhibitors agreed with the reported experimental data, and the 13b showed the most efficient binding to SARS-CoV-2 3CLpro. From rational drug design strategies based on 13b, the polar moieties (e.g., benzamide) and the bulky N-terminal protecting groups (e.g., thiazole) should be introduced to P1' and P4 sites in order to enhance H-bonds and hydrophobic interactions, respectively. We hope that the obtained structural and energetic information could be beneficial for developing novel SARS-CoV-2 3CLpro inhibitors with higher inhibitory potency to combat COVID-19.